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Abstract 
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models of strongly correlated particles with single-particle correlated and un- 
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Integrable strongly correlated electron systems (see, e.g. [1]) have been an important 
subject of research activity in recent years since they are believed to play a promising 
role in unraveling the mystery of high-Tc superconductivity. Several integrable correlated 
fermion systems which manifest superconducting properties so far appeared in the literature. 
Most famous are the supersymmetric t — J [2-4] model and the Hubbard model[5]. Other 
integrable correlated electron systems of interest include the correlated hopping model [6] 
and its generalization extensively investigated in [7,8]. The models considered in these last 
papers describe the dynamics of two type of particles (spin up and spin down) with kinetic 
terms given by correlated single particle hopping and uncorrelated hopping in the case of 
pair motion. 

In this letter, we introduce two new integrable quantum chains with correlated single- 
particle and uncorrelated pair hopping but having only one type of particles. 

We suppose that particles on the chain may be isolated (both nearest neighbour sites 
are empty) or part of a two-atom molecule (one of its nearest neighbour sites is occupied 
and the other one is empty) . This supposition is achieved by imposing the restriction of no 
simultaneous occupancies of three nearest neighbor sites. It means that two particles on the 
neighbor sites create the two-atom molecule which can hop as a whole or disintegrate. The 
Hamiltonian of such system in the most general form can be presented as follows 



H 



3 



X 



[(1 - nj-i){l - nj+2) + tinj-i{l - nj+2) 



+ 



t2{l - nj-i)nj+2 + ti2nj-inj+2] 



+ 



+ 



nj-inj+i[Vii + Vi2nj+2 + V2inj-2 



+ 



V22nj-2nj+2]}V 



(1) 



where crj~{aj ) creates (annihilates) a particle at site j and 



2 



is the corresponding occupation number. The model contains correlated single-particle hop- 
ping, which is described by the parameters ti, ^2, ^12 and uncorrelated pair hopping described 
by the parameter tp. We have also two- , three- and four-particle static interactions between 
nearest-neighbors. The operator V in (1) projects out any configuration with simultaneous 
occupancies of three nearest neighbor sites and we assume periodic boundary conditions. 

Certainly the Hamiltonian (1) is not integrable for an arbitrary choice of parameters. 
Therefore our problem is to find under what conditions the Hamiltonian (1) can be treated 
by the coordinate space Bethe ansatz technique. 

The obvious way to identify a state of the Hamiltonian (1) with n particles is to specify 
their positions Xi, . . . ,Xn ordered so that 

l<Xi<X2< ...<Xn<N. (2) 

We assume the following ansatz for the wavef unction. If we have only particles (xj+i 7^ 
Xi + 1, i — 1,2, . . . ,n — 1) we write 

n 

^{xi,...,Xn) ^Y. ^Pi'.^.Pn exp(« Yl ^Pj^j) (3) 
P 3=1 

where P is the permutation of 1, 2, . . . , n and kp are unknown quasiparticle momenta. The 
n superscript 1 in the amplitude indicates we have only isolated particles. In the case we 
have a pair at the position xi, — xi-\-l) the ansatz is 

'^{X'i_,...,Xi,Xi+i,...,Xn) = 
n 

E ^k:'.p;p,+...p. exp(z ^ kp^x,) (4) 
p j=i 

where the bar at the l^^ and (/-|- l)*'^ position of the superscript indicates the pair's position. 
The general case with many isolated particles and pairs follows from eqs.(3) and (4). 

It is not difficult to consider the eigenvalue equations if n = 1, 2 or for general n in 
the case we have only isolated particles. These equations give the connections between the 
coefficients A\,:^]p^. 
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where 



*PlP2 — ^^PlP2^PlP2 ^^PlP2^PlP2 

N^p^ = 1 - e*^^2l^ii + e*('=^i+*=^2) 

_l_ g«(2fcpi+fcp2) + g«(fcpi+2^P2) — tpe^**-^^i''''^^2). (6) 

In contrast to the XXZ model [2,9], this is not sufficient to prove that the Bethe ansatz 
works. In order to do that we must consider the eigenvalue equations at the boundary of the 
inequahties (2) for the case n — 3 and 4. This gives a complicated system of equations for 
the coupling constants of the Hamiltonian (1) . We have treated this system on a computer 
and found the following integrable cases 

(1) ti = ^2 = 2coshr;, ti2 = 1, tp = -e, 

T^i = 1/22 = 0, [/ = 2ip; l/i2 = 1^21 = e{ti - 1), (7) 

and 

.^s sinh27? ,„ sinh2r? £ cosh 77 

cosh 377 cosh 377 cosh 377 

Vu = 2tp cosh 27;; U = 2tp{l + 4 sinh^ r) cosh 27;) ; 

Vu = (e"^'' - 2 cosh 277)^^; V21 = (e^'' - 2 cosh 277)^^; 

V^22 = 0, (8) 

where e — ±1. 

The Bethe-ansatz equations are derived following the standard procedure [1-2]. Each 
state of the Hamiltonian is specified by a set of particle rapidities Xj{j — 1, . . . , 7i) related to 



the momenta of particles kj. The rapidities have to satisfy the Bethe-ansatz equations. For 
the model (7) these equations have the following form 



e sin(Aj + irj) 
sm{Xj — ifj) 



N-n 



^ g-ip -TT sin(Aj - A; + i7]) 
IJ^ sm{Xj - Xi- if)) 
sin(Aj + if]) 



sin(Aj — if)) 



(9) 
(10) 



and for the model (8) these equations are 



e sin(Aj — irj) 
sin(Aj + ir)) 



N-n 



_ _ip -A- cos(Aj - A; + i7]) sm{Xj - Xi - 2ir]) 
f-J^ cos(Aj — A; — irj) sin(Aj — A; + 2irj) 
sin(Aj — irj) 
sin(Aj + 17)) 



(11) 
(12) 



In both cases 



1=1 

is the momentum and the energy of the system is given by 

sinh^ 2?7 



(13) 



cosh 2r) — 



(14) 



cosh 2r] — cos 2Aj 

It is interesting to observe the similarity of the Bethe ansatz equations (9) and (11), and those 
of the spin-1 Zamolodchikov-Fateev model [10] and Izergin-Korepin model[ll], respectively. 
Although both models are exactly integrable , let us restrict ourselves to the more physically 
interesting model (7) with density p < 1/2. There are different regions, with distinct physical 
properties depending on the parameter A = —£ cosh 77. 

1) A =< —1, e = +1, 1] is real . We have a similar situation as in the ferromagnetic non- 
critical X X Z-chedn with fixed magnetization [2,12]. The ground state must contain exactly 
one string of maximum length. It means that there is a gap for the arbitrary concentration 
of particles and the system is in a phase separated region, where all the particles prefer to 
stay together. 

2) A>1,£=— l,77is real. We may analyse this case by considering the hmiting case 
7] — >■ +00. From this analysis it is clear that the ground state contains n/2 bound pairs 
characterized by a pair of complex particle rapidities 



A^ = ^(Ai^^±^r7). (15) 

Inserting (15) in eqs.(lO) and introducing the density function p(A) for the distribution of 
A^^-* in the thermodynamic hmit, we obtain the hnear integral equation 



27rp(A) + [2$(i) (A - A'; 2r]) + (A - A'; 47]MX')dX' 

— Ao 



ii-^m^'Hx,v)+^^'H\m, (16) 



where 



^"'(A,.)= , . (17) 

cosh 7] — COS A 

The parameter Aq is determined by the subsidiary condition for the total density p — n/N 
of particles 

j\{\)d\ = ]^p. (18) 
Similar calculations gives us the ground-state energy 

= 2p cosh 277-2 sinh 2r] t° [^^^^ (A, rj) + ^^^^ (A, 3ri)]p{X)dX. (19) 

iV J—Xo 

The other two regions are obtained from — 1 < A = — £ cosh 77 < 1. Denoting rj = i'-f,X — 
ip and instead of (9) we have 

£sinh(//j+i7) ^ ^_.p " sinh(/Xj - pi + i^) 
sinh(^j — i^) f-J^ smh.( pj — pi — ij) 

3) In the case — 1<A<0,0<£ = +1, 0, 7 < 7r/2 we solved numerically eq.(20) for 
up to 100 and checked that state now contains only strings of minimal length 1, [13] i.e., 
all particle rapidities {pj} have imaginary part i7r/2 (antiparticles). In the thermodynamic 
limit we then have the following integral equations for the distribution function 



27ra{p)- r ^'-'\p-p',2^)a{p')dp' 

J— no 



n 



{l-^)<^^'\p,2^), (21) 
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where 



a{i,')di,' = p, (22) 

-«n 



^'-'(^■7)- ■ ,-7" ■ (23) 

cosh jj, — cos 7 cosh jj, + cos 7 



and the ground-state energy is given by 

= -2p cos 27 + 2 sin 27 f^' ^^^^ (/x, 2-f)a{ii)di^. (24) 



4) In the case 0<A<1,£=— l,0<7<7r/2 our numerical results of (20) for lattice 
sizes N up to 100 indicate that the ground state contains only bound pairs, characterized 
by a pair of complex particle rapidities [14,15], like in region (15) 

=^(/^a±^7) (25) 
For the density function cr^^^ (//) we have the following integral equations 

2na^'\^i) + /';j2$(2)(/. - 27) + - ii'-Al)W'\li')dix' = 

(l-p)[$(2)(/.,7) + *(')(/^,37)] (26) 



J_\^^a^^\t,W-lp, ^^^^ 



(2) 

-E = 2pcos 27 - 2 sin 27 [$(2)(//, 7) + 37)]<7(^)(//)d/i. (28) 



(2' 

and the ground-state energy is given by 

(2) 
^0 

(2 

-Mo 

Solving numerically the corresponding integral equations in regions 2, 3 and 4 we show in 
Fig.l the ground-state energy as a function of density for some values of A. Except region I 
in all regions we expect gapless excitation for p < 1/2 . In order to understand the physical 
properties of the model under consideration we shall investigate the long-distance behavior 
of the correlation functions. For this purpose we shall use two-dimensional conformal field 
theory [16,17] and analytic methods [18] to extract finite-size corrections from the Bethe- 
ansatz equations. The results of these calculations indicate that the critical fluctuations 



are described by a conformal field theory with the central charge c — 1. The long-distance 
power-law behavior of the density-density correlation functions is given by the general form 

(p(r)p(O)) ~ p2 + Air-^ + Asr"" cos{2kFr)] 2kF = 7rp; (29) 



(30) 

while the pair correlation function is given by 

Gpi^) = (^/^Z+i' ^7+r^7+r+i) - Br-^. (31) 

The exponents a and (3 describing the algebraic decay are calculated from the dressed charge 
function = ^(''^o) 

This function is obtained by the solution of the integral equations (16), (21) and (26) with 
the right-hand side replaced by (1 — p). In Fig. 2 we show the exponent /3 as a function of p 
for several values of A in regions 2 and 4. Our numerical results indicate that as p — > 1/2 
the exponent j3 tends toward the value A'^/ti and /3 = 8(7r — 27)/7r in regions (2) and (4), 
respectively. In the region with dominant density-density correlations (5 > 1 the particles 
prefer to move individually, instead by pair hopping, but in the region with dominant pair 
correlations, /3 < 1 they create two-atom molecules. For arbitry values of A > 1 there 
exists a curve A = Ao(p) separating both behaviors. An analogous behavior of correlation 
functions can be observed in the models [7,8] which is translated in a strong tendency to the 
superconductivity. In region 3 we have no pairs . In Fig. 3 we show the exponent j3 which 
describes now the spin-spin correlation function < a^ajj^^ >. Our numerical results indicate 
that as p ^ 1/2 the exponent /3 tends toward the value P — A'j/n. 
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FIGURES 

FIG. 1. The ground-state energy as a function of the density p for some values of A = — £ cosh ry, 
in region 2,3 and 4. a) A = 11.5920, b)A = 2, 5092, c)A = 0.7071, d) A = 0.1423, e) A = -0.9239, 
f) A = -0,8090, g)A = -0.5 

FIG. 2. The exponent (3 describing the pair -pair correlation function as a function of the 
density p for some values of A = —e cosh 77 , in regions 2 and 4. a) A = 11.5919, b)A = 1.2039, c) 
A = 0.7071, A = 0.3827, e) A = 0.1423, f)A = 0.0383 

FIG. 3. The exponent /9 describing the spin-spin correlation function as a function of the density 
p for some values of A = — ecoshr/ , in region 3. a) A = —0.5, b) A = —0.7071, c) A = —0.8090, 
A = -0.9239, e)A = -0.9511, f)A = 0.9980 



12 



